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Abstract 

I: 

equations of pions and nucleons in low temperatures and low baryon number densities. The n to 
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entropy density ratio n/s maps out the nuclear gas-liquid phase transition by forming a valley 
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tracing the phase transition line in the temperature-chemical potential plane. When the phase 
transition turns into a crossover, the n/s valley gradually disappears. We suspect the general 
' feature for a first-order phase transition is that n/s has a discontinuity in the bottom of the n/s 

o ' 

CO ■ valley. The discontinuity coincides with the phase transition line and ends at the critical point. 

Beyond the critical point, a smooth n/s valley is seen. However, the valley could disappear 
further away from the critical point. The n/s measurements might provide an alternative to 
' identify the critical points. 
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I. INTRODUCTION 



Shear viscosity rj is a transport coefficient which has recently been attracting lots of 
attention. It characterizes how strongly particles interact and move collectively in a many- 
body system. In general, the stronger the interparticle interaction, the smaller the shear 
viscosity. It is conjectured l| that, no matter how strong the interparticle interaction is, 
the shear viscosity to entropy density s ratio has a minimum bound 1 /An. I.e. rj/s > 
in any system. The bound was motivated by the uncertainty principle and the observation 
that n/s = 1/47T for a large class of strongly interacting quantum field theories whose dual 
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descriptions in string theory involve black holes in anti-de Sitter space [2l,l3|,y,l5|. In Ref. 

supporting evidence of the conjecture was given for matters like H 2 0, He and N. Their 
rj/s curves reach their minima near the gas-liquid phase transitions with the bound well 
satisfied. Recently, rils close to the minimum bound was found in relativistic heavy ion 



collisions (RHIC) 



0,0, Is 



(and in lattice simulations of a gluon plasma 



9|) just above the 



deconfinement temperature T c (~ 170 MeV at zero baryon density [lOj). This suggests 
the quark gluon plasma (QGP) is strongly interacting at this temperature, which is quite 
different from the traditional picture of weakly interacting QGP. * (However, see Ref. 
for a different interpretation.) Also, rj/s close to the minimum bound was found in cold 
fermionic atoms in the infinite scattering length limit 20J. A relation between rj/s and 



the jet quenching parameter in QGP was proposed in Ref. 
In Refs. 



221 and 



2l| 



231 ]. it was found that rj/s of QCD in the confinement and decon- 



finement phases is qualitatively different. When T < T c (the confinement phase), rj/s 
is monotonically decreasing in T because the system is dominated by goldstone bosons 
which interact more weakly at lower T. When T > T c (the deconfinement phase), rj/s 
is monotonically increasing in T because the interaction between quarks and gluons is 
weaker at higher T due to asymptotic freedom. It makes perfect sense to have the phase 
transition from the point of view of preserving the 77/s minimum bound. This is because 
if the qualitative behavior of 77/5 is not changed by a phase transition (or crossover), then 
the bound could be violated. One concludes that the minimum or valley of the r)/s curve 
lies in the vicinity of T c (the extrapolation of the low(high) temperature rj/s curve sets a 
upper (lower) bound on T c ) 23]. This behavior is also seen in the H 2 0, He and N systems. 



* See als o Ell . H3 , fl^ ] . For discussions of the possible microscopic stucture of such a state, see 0, 0, 



16, 



13,1131 



It was further noticed that below the critical pressure, a cusp appears at the minimum of 
r]/s which coincides with the critical temperatures 



22 . 



In Ref. 24|, it is argued that a universal minimum bound on rj/s should not exist. The 
counterexample given is a system of mesons made by heavy quarks and light antiquarks. 
Since s scales linearly with the number of heavy quark flavor Nf and rj is insensitive to 
Nf, the rj/s > 1/Att bound could be violated in some special large Nf limit. However, 
this system is metastable. As far as the qualitative relation between T c and the valley of 
rj/s is concerned, it does not matter the minimum bound of rj/s is 1/Att or 0. As long as 
there is a lower bound, the monotonic behavior of r]/s will be affected by the bound. 

In this manuscript, we extend the discussion of the r)/s of QCD in the confinement 



phase at zero baryon chemical potential /i 23( to finite [i and study its relation to the 



QCD phase diagram. String theory methods give f]/s = 1/47T for M = 4 supersymmetric 
theories with finite R-charge density, suggesting the minimum bound is independent of 
/i. For QCD, the fermion sign problem (the fermion determinant is not positive definite) 
makes the current lattice QCD methods inapplicable in the low T and finite [i regime. An 
alternative is to use effective field theory (EFT). Reliable results using EFT in hadronic 
degrees of freedom can be obtained when both T and /i are small. At higher /i (with 
\kpa\ ^> 1, where kp denotes the fermi momentum and a is the nucleon-nucleon scattering 
length) the problem becomes non-perturbative in coupling and mean-field treatments are 
not sufficient. (It is essentially the same type of problem as in cold fermionic atoms 
near the infinite scattering length limit.) Non-perturbative computations of the EFT on 
the lattice is free from the fermion sig n problem at the leading order (LO) with only 
non- derivative contact interactions 25[. But this theory is suitable only in low T and 
low density systems. For the nuclear matter problem, the inclusion of one pion exchange 
will re-introduce the sign problem. However, the sign problem is claimed to be mild and 



lattice simulations are still possible 26] . The computation of r] using lattice nuclear EFT 
has not been carried out before. Although 77 is associated with real time response to 
perturbations, it can be reconstructed through the spectral function computed on the 
Euclidean lattice @, Q. 

As an exploratory work, we compute 7] using coupled Boltzmann equations for a system 
of pion 7r and nucleon N while the entropy s is computed only for free particles. This 
approach will not give accurate rj/s in the regime dominated by near threshold NN 
interaction. However, for most of the regime we are exploring, our result should be 



T, GeV 



QGP 




FIG. 1: A semi-quantitative sketch of the QCD phase diagram [28[ courtesy of M. Stephanov. 
robust. 



II. LINEARIZED BOLTZMANN EQUATION FOR LOW ENERGY QCD 

We are interested in the hadronic phase of QCD with non-zero baryon-number chemical 



potential \i. For reference, a schematic QCD phase diagram from a recent review 28j] is 
shown in Fig. 1. The detail structures of the quark matter phases are still unclear. 
When /i<mjv (= 938 MeV, the nucleon mass) the nucleon population is exponentially 
suppressed. The dominant degrees of freedom are the lightest hadrons — the pions. The 
pion mass m 7T (= 139 MeV) is much lighter than the mass of the next lightest hadron — 
the kaon whose mass is 495 MeV. Given that T c is only < 170 MeV, it is sufficient to 
just consider the pions in the calculation of thermodynamical quantities and transport 
coefficients for T <C T c . When — // < m n , the nucleon population is no longer 
suppressed compared with the pion. We will limit ourselves to the low T and low \i region 
where only ir and N are important degrees of freedom. 

The shear viscosity of a system is defined by the Kubo formula 

V = -\ f df f dt I d* 3 ([T^(0),T^(x,t)]> , (1) 

^ J —OO J —OO J 

with T*- 7 the spacial part of the off-diagonal energy momentum tensor. The Kubo formula 
involves an infinite number of diagrams at the LO even in the weak coupling 4 theory {29]]. 
However, it is proven that the summation of LO diagrams in a weak coupling 4 theory 
is equivalent to solving the linearized Boltzmann equation with temperature dependent 



particle masses and scattering amplitudes |29j. We will assume the equivalence between 
the Kubo formula and the Boltzmann equation still hold in our ttN system. Later we 
will check whether the mean free path is still much larger than the range of interaction. 
This is a requirement to apply the Boltzmann equation which makes use of semi-classical 
descriptions of particles with definite position, energy and momentum except during brief 
collisions. 

In the Boltzmann equation of our irN system, the evolution of the isospin averaged 7r 
and N distribution functions f n ' N = f n ' N (x., p, t) = fp ,N {x) (functions of space, time and 
momentum) are caused by interparticle (tttt, ttN and iViV) collisions 
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+9N / ^{/^JiJ-Wfl , 

J 123 

^ffix) = fJ w dT^{fn?F?F?-F»F»f»f?} 
+9* [ dT% p {f?f?F?F? - F?F 2 N f*f p N } , 

J 123 



(2) 



(3) 



where F^^ = 1 ± /, 
factors g n = 3 and g N = 4. 
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where 7^v is the 7riV scattering amplitude with momenta 1,2 — > 3, p. The 7T7t and NN 
weighted integration measures are given analogously. We use the LO chiral perturbation 
theory (xPT) |30|, |3l|, |32j, |33j result for the isospin averaged im scattering amplitude in 



terms of Mandelstam variables (s, t, and u) 

1 



l^l 2 = 7777 [ 21 ™t + 9s 2 - 2Amls + 3(t - w) 2 ] , 



(5) 



where /„. = 93 MeV. The isospin averaged iViV scattering amplitude are described by 
effective range expansion 34j. In the center of mass (CM) frame 



\T NN \ 2 = 3(47rmjv) 5 



1 1 2 • , 

h ~r\P - ip\ 
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+ h 7T r 3P - «P 

a 3 2 
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where p is the magnitude of nucleon momentum in the CM frame and ai(3) and 7*1(3) are 
the spin singlet (triplet) scattering length and effective range, respectively, ai = —17.9 



a>3 = 5.42, r\ = 2.77, and r% = 1.76, all in units of fm. Note that near the threshold 
(p = 0), \T NN \ 2 is proportional to a\ + a\ which is greatly enhanced by the large scattering 
lengths. The interaction is smaller away from the threshold. 

The 7iN scattering [7r(gi) N(pi) — > 7r(g 2 ) N(p 2 ) } amplitude is also given by the LO 
XPT |35j. In the CM frame 

\% N \ 2 = ^(E N + m N ) 2 [2\g_\ 2 + q 2 sin 2 6\h + \ 2 ] , 

where gA = 1-26 is the nN coupling constant, En is the nucleon energy, q is the magnitude 
of pion momentum, 9 is the angle between qi and q 2 , and u is the pion energy. The thermal 
corrections for ttti, ttN scattering amplitudes and particle masses are higher order in %PT. 
The T dependence in NN scattering amplitude is also small at low T. Because the thermal 
correction for the inverse scattering length scales as T which is much smaller than the 



thermal momentum ~ ^/itlnT^ it can be neglected. 

In local thermal equilibrium, the distribution functions are fl{x) = (e^ x ^ v '* t ^ pli — l) 
and J p (x) = ( e /%)WP + iy\ w here (3{x) is the inverse temperature, V^(x) is the four 
velocity of the fluid at the space-time point x and p M = (£r — /i, p) in the V(x) = 
frame. A small deviation of f p from local equilibrium can be parametrized as 

f p (x)=f p (x)[l-F l p (x)4(x)] , 1 = tt,N, (8) 
with F; v v = 1±/: v '. The energy momentum tensor is 



9wfp(x) 9Nfp(x) 

E 77 E N 
p p 



(9) 



We will choose the frame with zero fluid velocity ~V(x) = at the point x. This implies 
d v V° = after taking a derivative on V lJL {x)V^{x) = 1. Furthermore, the conservation 
law at equilibrium d^T^ u \ Xv= Q = allows us to replace d t (3{x) and <9 t V(x) by terms 
proportional to V • V(x) and V/3(x). Thus, to the first order in a derivative expansion, 



t We have used the isospin averaged scattering length defined as a\ = (a^„ + a 1 + a? , ) /3 



X p ( x ) can be parametrized as 



= A'(p)V • V(x) + ^.(p) - - -^V • V(x) J , (10) 

= B\p) (pm - U l3 

where i and j are spatial indexes. A and 5 are functions of x and p. However, we have 
suppressed the x dependence. 

Substituting Eq. fflOl) into the Boltzmann equation Eq. (T5]), one obtains a linearized 
equation for B 71 : 

(p iPj - ~w) = ^ / 23 

[B«.(p) + B*(k 3 ) - B^h) - B^h)] 

J 123 

x [2$(p) - Bf^h) + B£ (A*) - B» (hi)] 
= g^G™ [B-} + g N G x f m + g N G 2 f [B N ] . (11) 

The analogous equation for B N is obtained by replacing tt <-> A" on the right-hand side 
of the above equation. There are two other integral equations involving A l (p)V ■ V(x) 
which are related to the bulk viscosity (. They will not be treated in this work. 

In fluid dynamics the energy momentum tensor at equilibrium depends on pressure 
P(x) and energy density e( T$\x) = {P(x) + e(x)} Vp{x)V v {x) - P{x)8„ v . A small 
deviation away from equilibrium gives additional contribution to T^ u = T$ + 5T^ U whose 
spatial components define the shear and bulk viscosity 



x 



(12) 



Comparing the above definition with Eqs. ([9]) and ( [101) . we obtain 



r, = g,L'\B'}+ 3N L N [B N ] , 

Now one sees immediately that if all the tttt, tcN and A^A" scattering cross sections are 
reduced by a factor A, then Eq. ffTTj) implies the B l functions will be A times larger. Then 
by Eq. ffl3|) . r) will be A times larger as well. This is the non-perturbative feature of the 
Boltzmann equation. It gives a divergent r\ for a non-interacting theory. 



Contracting both sides of Eq. fflTl) by (pipj — %/3) and applying it to Eq. ffT3l) yields 

1 10 J (2TifE; Jp p l3KP > 

x {gJJ% [B"] + gsGrf [B"\ + g N G 2 f [B N ] } 

+ [tt «-> N] 

= gl{B n \G™ [B v ]) + g v g N {(B v \G x * N [B v ]) + (B^\G^ N [B N ])} 

+ [tt <-> iV] . (14) 

To compute 77, one can just solve B™( N \p) from Eq. (fTTl) . But here we follow the 
approach outlined in Ref . 36|, [37| to assume that B n ^ (p) is a smooth function which 
can be expanded using a specific set of orthogonal polynomials: 



B\p) = |p| tf f> r S/ P) (*(p)) 

r=0 

00 

= J2^\z(p)) , I — 7T, N 



(15) 



r=0 



where B^ r \z) is a polynomial up to z r and 6[ is its coefficient. The overall factor |p|^ will 
be chosen by trial and error to get the fastest convergence. We find that using y = 1.89 
and z(p) = |p|, the series converges rather rapidly. The orthogonality condition 



£ f d 3 p |p| 2+ ^ ^ 
15 J (2vr) 3 E l Tp r 



T p t p b!t\z)bM(z) = l?>&. 



(r). 



(16) 



can be used to construct the B\ r \z) polynomials up to normalization constants. For 
simplicity, we will choose 

(17) 



B?\z) = l. 



With this setup, Eqs. flI3j) and §Fh yield 



£ LKL^ + g N b?L% 



} 0,r , 



while Eq.flHD yields 



J2 9*K {a^s (bP\GT \b^] ) + g N fe (b^\G^ n \B^ 



r,s=0 



+ [tt «-> N] . 



b? B^\G 2 * N 



B 



iV 
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The b^s and are functions of li and T. In general they are independent of each 
other. Thus the above two equations yield 

where the matrix elements are given by 



/ . \ A^vr rs NN rs 



n7l rs = (bW\g™[bP]) + ^(bP\gz n [bP]) , 

nN rs = (BP\G% N [B$]) , 

NN rs = (BP\G NN [Btf}) + ^(BP\G?*[B$}) , 

9n 

N7t rs = (bP\G%*[BW]) . (21) 

The matrix equation in Eq. (1201) allows us to solve for b s and obtain the shear viscosity. 
Since the expansion in Eq. (fl5|) converges rapidly, one does not need to keep many terms 
to solve for rj. If only the s = term is kept, then 

_l / r(0) \ 

/ Lilt 1 

:22i 

J N 

The computation of the entropy density s is more straightforward since s, unlike rj, does 
not diverge in a free theory. The contributions from hit and 7rN scattering are perturbative 
and are higher-order effects in %PT. The NN contribution is also perturbative at low li, 
but it becomes non-perturbative at li « itln and low T. In this regime, the system is 
governed by near threshold NN interaction. With that limitation in mind, we compute 
s as a free gas: 

R 2 d <h z -k + 9n log Z N 
S ~ ~^ dp ' {26) 

where the partition functions Z^n) for free pions(nucleons) is 

logZ 7r(JV) _ 1 f d 3 p ; i _ 
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v - 1 







TTTT 00 TTiV 00 \ / Li 0) 

iVvr 00 NN 00 j \ L$ 



log lTe"^ , (24) 
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with Ep = E£ and E N = E N — li, up to temperature independent terms. 



III. rj/s AND THE QCD PHASE DIAGRAM 



In Fig. 2, rj/s as a function of T and /i is shown as a 3-D plot and a contour plot. Note 
that at the corner of large \x and large T, the system is no longer in the hadronic phase. 
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FIG. 2: r//s of QCD shown as a function of T and fi shown as a 3-D plot (a) and a contour 
plot (b). Note that the corner of large fi and large T should be discarded since it is not in the 
hadronic phase. 

Thus, the result should be discarded there. In general rj/s is decreasing in T except when 
fi ~ mjy and T < 30 MeV. (This regime is blown up in Fig. 3 and will be studied later.) 
There are some interesting structures at larger fi, but rj/s is fi independent when fi < 500 
MeV . This is because when fi <C mjy the nucleons only exist through particle-antiparticle 
pair creations, thus they are highly suppressed. The rj/s is determined by the pion gas 
which isu independent. Our result just reproduces the fi = result of Ref. [23J (see 



36, 



37 



38 



39( for earlier results) in this regime. 
When fi > rriN — m n = 800 MeV, the nucleon population is no longer suppressed 
compared with the pion population and when fi > mjy the nucleons become the dominant 
degrees of freedom. Numerically rj/s is dominated by the nucleon contributions when 
fi > 800 MeV. It is decreasing in both T and fi until fi ~ m^. This is because s 
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FIG. 3: rj/s of QCD near the nuclear gas-liquid phase transition shown as 3-D plots (a) and 
(b), viewed from different angles, and a contour plot (c). The rj/s maps out the nuclear gas- 
liquid phase transition shown in Fig.l by forming a valley tracing the nuclear gas- liquid phase 
transition line in the T-ji plane. When the phase transition turns into a crossover at larger T, 
the valley also gradually disappears at around 30 MeV. There should be a discontinuity that 
looks like a fault in the bottom of the rj/s valley that is not seen in our approximation. The fault 
would lie on top of the phase transition line and end at the critical point (T ~ 10-15 MeV in our 
result). Beyond the critical point, rj/s turns into a smooth valley. The valley could disappear 
far away from the critical point. Similar behavior is also seen in water shown in Fig. 4. We 
suspect these are general features for first-order phase transitions. 

is increasing in both T and ii while rj is getting smaller at higher ii (larger nucleon 
population) and lower T (stronger interaction, closer to the interaction threshold). 500- 
800 MeV in /i is the transition between the tt and iV dominant regimes. 

Now let us focus on the /i ~ and T < 30 MeV region in the rj/s plot. Two 3-D plots 



viewed from different angles are shown in Fig. 3(a) and 3(b) and a contour plot is shown 
in Fig. 3(c). One clearly sees that rj/s maps out the nuclear gas-liquid phase transition 
shown in Fig.l by forming a valley tracing the nuclear gas-liquid phase transition line in 
the T-fi plane. When the phase transition turns into a crossover at larger T, the valley 
also gradually disappears at around 30 MeV. This result is encouraging. 

However, even though the gross features of the phase transition are mapped out by the 
rj/s valley nicely, some details are not correct. First, since the density is discontinuous 
across the first-order phase transition, rj, s, and rj/s are likely to be discontinues across the 
phase transition as observed in H2O, He and N systems in Ref. [221 ] . This discontinuity, 
which defines the critical chemical potential ji c , should lie at the bottom of the rj/s valley. 
Second, the position of fi c suggested by our result is not correct. Near T = 0, one expects 
ri c ~ m N — (B), where (B) is the binding energy per nucleon, but we have ji c > m N . 

It is quite obviously that our free particle treatment of s is very poor near the phase 
transition. However, we have not pursued other treatments like the mean field approx- 
imation in this work because it is known that the approximation is insufficient when 
|/cfo| ^> 1- For the same reason, the computation of 1] using the Boltzmann equation 
might not be justified near the bottom the valley even though the mean free path is 
still bigger than the range of potential (~ 1 fm). However, the valley of 1] is located at 
ji < m N near T = 0; thus, it is possible that after reliable s is used, ji c for rj/s will 
be in the correct position. Furthermore, the regime in Fig. 3 is completely dominated 
by the nucleon degree of freedom (rj/s hardly changes with the thermal pions completely 
ignored) . This simplifies the problem significantly and makes the system exhibit universal 
properties shared by dilute fermionic systems with large scattering lengths such as cold 
atoms tuned to be near a Feshbach resonance. 

As mentioned above, it was observed that below the critical pressure, rj/s has a dis- 
continuity at the critical temperature for H 2 0, He and N, and above the critical pressure, 
rj/s has a smooth minimum near the crossover temperature (defined as the temperature 



where the density changes rapidly) [22j . For QCD at pi — 0, rj/s also has a valley near the 
crossover temperature 22|, |23J. But there is no evidence yet to show the valley is smooth 
or has a discontinuity. 

We suspect rj/s should be smooth in a crossover and should have a discontinuity across 
a first-order phase transition. If this is correct, then in nuclear gas-liquid transition, there 
should be a discontinuity looking like a fault in the bottom of the rj/s valley. The fault 
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FIG. 4: 77/s of water shown as a function of temperature and pressure. Below the critical 
pressure 22.06 MPa, 77/s has a discontinuity at the bottom of the valley. Above the critical 
pressure, the valley becomes smooth and the bottom (minimum) of the valley moves toward 
the large T and large P direction. In a limited range of T, 77/s could look like a monotonic 
function without a valley far away from the critical pressure. The similar behavior is seen in all 
the materials with data available in 



43 



44J, including Ar, CO, C0 2 , H, He, H 2 S, Kr, N, NH 3 , 



Ne, O, and Xe. 



would lie on top of the phase transition line and end at the critical point where the 
first-order phase transition turns into a crossover. If we look at Fig. 3(c), we would 



conclude that the critica 
experimental extractions 



point is at T ~ 10-15 MeV, agreeing with 7-16 MeV from 
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41 



42j. Near the critical point, a smooth 77/s valley is seen 



in the crossover (like the confinement- deconfinement crossover of QCD at // = 0); however, 
the valley could disappear far away from the critical point. We suspect these are general 
features of first-order phase transitions. Indeed, this behavior is seen in all the materials 



with data available in the NIST and CODATA websites [43j,|44j, including Ar, CO, C0 2 , 
H, He, H2O, H2S, Kr, N, NH3, Ne, O, and Xe. As an example, we plot the 77/s of H2O 
as a function of T and P in Fig. 4. Below the critical pressure 22.06 MPa, 77/s has a 
discontinuity at the bottom of the valley. Above the critical pressure, the valley becomes 
smooth and the bottom (minimum) of the valley moves toward the large T and large P 



direction. In a limited range of T, rj/s could look like a monotonic function without a 
valley far above the critical pressure. Thus, one might use rj/s measurements to identify 



22 



4q for a similar point 



the first-order phase transition and the critical point (see Refs. 
of view) . 

Also, for all the materials listed above except CO, the rj/s has a positive jump across 
the phase transition line from the gas to the liquid phase as shown in Fig. 4. For CO, the 
jump is positive at smaller pressure but becomes small and negative at higher pressure. 
It is interesting to recheck wether CO really is an anomaly. However, even without 
considering CO, the sign of the rj/s jump for first order phase transitions might not be 
universal, either. For QCD in the limit of a large number of colors, the jump is negative 
from the low to high temperature phases 



IV. CONCLUSION 



We have computed the shear viscosity of QCD in the hadronic phase by the coupled 
Boltzmann equations of pions and nucleons in low temperatures and low baryon number 
densities. The ratio rj/s maps out the nuclear gas- liquid phase transition by forming a 
valley tracing the phase transition line in the temperature-chemical potential plane. When 
the phase transition turns into a crossover, the rj/s valley also gradually disappears. We 
suspect the general feature for a first-order phase transition is that rj/s has a discontinuity 
in the bottom of the rj/s valley. The discontinuity coincides with the phase transition line 
and ends at the critical point. Beyond the critical point, a smooth rj/s valley is seen on 
the crossover side. However, the valley could disappear far away from the critical point. 
The rj/s measurements might provide an alternative to identify the critical points. 
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